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Based on the previous approach, we have investigated a possibility to resolve the discrep¬ 
ancy between the E710, E811 and CDF at ^/s = l.STeV, using the experimental data of the 
pp, pp total cross sections and ratio up to the SPS experiments (-^5 = 0.9TeV) as 
inputs. We predict and p^^ at the Tevatron energy (yT = l.STeV) as crf^t = 75.9±1.0mb, 
pPP = 0.136 ±0.005. 

It turns out that only the data of E710 is consistent with the prediction in the one 
standard deviation. So we can conclude that E710 is preferable but we can exclude neither 
CDF nor E811 results. 


§1. Introduction 

Recently,^) we have searched for the simultaneous best fit of the average of pp, pp 
total cross sections ( ), and the ratio of the real to imaginary part of the forward 

scattering amplitude( ) for 70GeV < Piab < Piarge as inputs in terms of high- 
energy parameters cq, ci, C 2 and ftp' constrained by the FESR with Af(~ lOGeV). 
Block and Halzen^)’^^ also reached to the similar conclusions independently based on 
duality in a different approach. We first chose Piarge = 2100GeV corresponding to 
the ISR region( ^/s ~ GOGeV ). Secondly we chose Piarge = 2x lO^GeV corresponding 
to the Tevatron energy( ^/s ~ 2TeV ). We then predicted and at the LHC 
and the high-energy cosmic-ray energy regions. It turned out that the prediction of 
agrees with pp experimental data at the cosmic-ray regions'^^^®^ within errors in 
the first case( ISR ). It has to be noted that the energy range of predicted p^^^ 
is several orders of magnitude larger than the energy regions of p^^^ input. If 
we use data up to Tevatron( the second case ), the situation has been much improved 
although there are some systematic uncertainty coming from discrepancy of the data 
between E7I0,^) ESll^^ and CDF®) at ^/s = l.STeV.^) Finally we concluded that 
the precise measurements of in the coming LHC experiments will resolve this 
discrepancy at ^/s = l.STeV. 

The purpose of this paper is to investigate a possibility to resolve this discrepancy 
using the experimental data of and up to the SPS experiments (-^/i = 
0.9TeV). 
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§2. The general approach 


As in the previous paper/^ let us first consider the crossing-even forward scat¬ 
tering amplitude defined by 


^with Jnt ^ "“A) 

We also assume 

Im = Im -R(z^) -|- Im Fpi{v) 

(cO + Cllog^ + C2log2^)+^( 


( 2 - 1 ) 


V 


V \ °‘p' 

M. 


(2-2) 


at high energies {u > N). It is to be noted that cq, ci, C2 and fdp/ are dimensionless. 
We have defined the functions R^n) and Fpi{v) by replacing /r by M in Eq. (3) 
of ref. 10). Here, M is the proton( anti-proton) mass and v, k are the incident 
proton(anti-proton) energy, momentum in the laboratory system, respectively. 

Since the amplitude is crossing-even, we have 


Ri^) = 


IV 

2M2 


2co -I- C27r^ -I- Cl log 


-h log— 

M 


+ C2(log^^+log^^) 


Fp'{v) = -■ 

and subsequently obtain 


I3p, ({e-'^^v/M^p' +{v/M)°‘p' 


M 


SinTTOp/ 


Re Fp.{„) =. 

^ ^ M \MJ 

substituting op/ = ^ in Eq. (( 231 ). 

EESR : The FESR corresponding to n = jg. 
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pN 
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Jo 


vim Fpi{v)dv 


(2-7) 


We call Eq. (12-71) as the FESR which we use in our analysis. 

The ratio; The ratio, the ratio of the real to imaginary part of F^F 

was obtained from Eqs. ( 1231 ), 1131 ) and ( 1231 ) as 

= 


fp-i. Re F^~^\v) Re R{v) + Re Fpi (v) 

' 17/1= - 


Im F(+)(i/) Im R{v) + Im Fpi{v) 
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55, (c.+2c2loga-%^ (*)“'" 

■ 

47r 

Although the numerator of Eq. (ESI) becomes large for large values of a real 
constant has to be introduced in principle since the dispersion relation for Re 
requires a single subtraction constant So, we also add F^^\0) in the 

numerator as 


3(+)(jy) = 


TTiy 

2M^ 


(ci+2Qlog^)-fe(j^)"» + F(+)(0) 

in 


(2-9) 


As will be discussed in the Appendix, the introduction of this constant slightly 
modifies the value of {v) although it will not affect the value of . So, we use 
the Eq. H2-9I) as the value of {v) in this analysis. 

The FESR, Eq. (l2-7j) . has some problem, i.e., there are the so-called unphysical 
regions coming from boson poles below the pp threshold. So, the contributions from 
unphysical regions of the first term of the right-hand side of Eq. (EiZI) have to be 
calculated. These contributions can be estimated to be an order of 0.1% compared 
with the second term.^^ Thus, it can easily be neglected. 

Therefore, the FESR, the formula of cr|^^(Eqs. (12-11) and (12-21) 1 and the p^^^ ratio 
(Eq. (1^ 1 are our starting points. Armed with the FESR, we express high-energy 
parameters cq, ci, C 2 ,/3p' in terms of the integral of total cross sections up to N. 
Using this FESR as a constraint for fipi = (5pi{co, ci, C 2 ), there are four independent 
parameters including E%^(0). We then search for the simultaneous best fit to the 
data points of (k) and p^^\k) for 70GeV< k < Piarge corresponding to the SPS 
energy {Piarge — 0.43 x 10®GeV {^/s = 0.9TeV)), to determine the values of cq, ci, C 2 
and giving the least x^- We thus predict the utot and p^^'^ in the Tevatron 

energy region (y^ = l.STeV). 


§3. Predictions for and at y^ = l.STeV 


Using the data up to y^ = 0.9TeV ( SPS ), we predict and p^~^^ at the 
Tevatron energy ( ^/s = l.STeV ). 

Analysis 1: As was explained in the general approach (§2), both and 
Re data in 70GeV < k < Piarge = 4:.3 x 10®GeV(y% = 0.9TeV ) are fitted 
simultaneously through the formula of ( Eqs. (12-11) and (12-21) 1 and the p^~^^ ratio 
( Eq. (12-91) 1 with the FESR ( Eq. (12-71) 1 as a constraint. 

The {k) data points are obtained by averaging and data points^^^ 
when they are listed at the same value of k. For the details of data treatment of 
and Re F^^\ see ref. 1). The FESR gives us 


8.87 = Co -|- 2.04ci -|- 4.26c2 -|- 0.367/3p' 


(3-1) 
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(Eq. (12) of ref. 1)), where we use the central value of ^ tot (k) = 3403 ± 

20GeV*) for N = lOGeV in Eq. (t^ . 

The result of the fit is shown in Fig. ^ The values of parameters and resulting 
are given in Tables HI and HIl respectively. 

Table I. The values of parameters in the best fit to the data up to SPS energy (y^ = 0.9TeV) in 
the analysis l(fit to the data in 70GeV < k < Piarge = 4.3 x 10®GeV ). The error estimations 
are done as follows: The C 2 is fixed with a value deviated a little from the best-fit value, and 
then the is done by three parameters co, ci and where dp' is represented by the 

other parameters through FESR(Eq. ll3Tll i. When the resulting is larger than the least 
of the four-parameter fit by one, the corresponding value of C 2 gives one standard deviation. 
The higher and lower dot-dashed lines in Fig. Q represent this deviation of C 2 . The errors of the 
other parameters are estimated through similar procedures. 



C2 

Cl 

Co /3p' 


Analysis 1 

0.0466 ± 0.0047 

-0.161 T 0.078 

6.27 ±0.33 7.45 T 0.51 

12.65 ± 5.66 


Table II. The values of x^ for the fit to data in 70GeV < k < Piarge = 4.3 x 10®GeV(Analysis 
1): Nf and Na{Np) are the degree of freedom and the number of data points in the 

fitted energy region. 



X^Nf 

Xl/Na xl/Np 

Analysis 1 

8.1/20 

5.7/17 2.4/8 


In terms of the best-fit values of parameters in Tabled the predictions at \/s = 
l.STeV are obtained as 

= 75.9 ± 1.0 mb, = 0.136 ± 0.005 , (3-2) 

where the errors correspond to the one standard deviation of C 2 , since the C 2 log^(z^/M)- 
term in Eq. (EH is most relevant for predicting in high energy region. (See the 
caption in Tabled) 

The equation EH has to be compared with the experimental values at y/s = 
l.STeV; 


(Tf^t(E811) = 71.71 ± 2.02 mb, 
aZ(B710) = 72.8 ± 3.1 mb, 

aZ(CnF) = 80.03 ± 2.24 mb, (3-3) 

where we note that the difference between and is negligible at the relevant 
energy. It is worthwhile to notice that only the data of E710^^ is consistent with the 
prediction, Eq. EH in the one standard deviation ( 72.8 -|- 3.1 = 75.9 ). 

If one tolerates two standard deviations, both CDF®^ (80.03 — 2.24 x 2 = 75.55) 
and E811®)(71.71 -F 2.02 x 2 = 75.75) are consistent with the predictions Eq. (13-21) . 

This value is obtained by numerically integrating the experimental k'^cy['^^ — +0'tot)/2. 

See, ref. 1) for details. 
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mb 



(c) : High energy region (d) 

Fig. 1. Predictions for and in terms of the Analysis 1. The fit is done for the data up 
to SPS energy, in the region 70GeV< k < 4.3 x 10®GeV(11.5GeV < < 0.9TeV) which is 

shown by the arrow. Total cross section in (a) all energy region, versus logioPia6/GeV, (b) 
low energy region (up to ISR energy), versus Piab/GeV and (c) high energy (Tevatron-collider, 
LHC and cosmic-ray energy) region, versus center of mass energy Ecm in TeV unit, (d) gives 
the Re /Im P^^^) in high energy region, versus Ecm in terms of TeV. The thin 

dot-dashed lines represent the one standard deviation of C2.(See the caption in Table ||]) The 
corresponding values of parameters are {c 2 ,c\,cq, j3pi ^ = (0.0466 ± 0.0047,-0.161 =F 

0.077, 6.27 ± 0.31, 7.45 T 0.48,12.65 ± 0.69). 
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So we can conclude that E710 is preferable but we can exclude neither CDF nor 
E811 results. 

The predictions at LHC energy (yT = 14TeV) in terms of the best fit values of 
high-energy parameters in Tableware 

o-fP = 107.2 ± 2.8 mb, = 0.128 ± 0.005 , (3-4) 

where the errors correspond to one standard deviation of C 2 . We should note that 
Eq. H3-4|) is consistent with the recent prediction by Block and Halzen,^) = 
107.3 ±1.2 mb, = 0.132 ± 0.001. 

An interesting observation: We can make the following interesting observation. 

We fitted the data for and above 70GeV, as is shown by the arrow in the 
Fig. 1(a), Fig. 1(d) to predict higher-energy data. It is interesting to observe that the 
prediction of are also in good agreement with experiments, even below 70GeV. 
The reason is as follows: The requirement of FESR, Eq. (EZI) is nearly equal to 
require that the theoretical value of is nearly equal to the experimental value 
at the upper limit of the integral N = lOGeV since higher side of the integral is 
enhanced because of k'^ in the integral. 

Because of this observation, we can apply the same formula to fit the data in 
the lower energy region than in the analysis 1. 

Analysis 2: Data in lOGeV < k < Piarge = 4:.3 x 10^GeV(4.54GeV < y/s < 
0.9TeV ) are fitted through the same formula in the analysis 1. Additionally 15(2) 
data points are included in {Re 

The result of the fit is shown in Fig. The values of parameters and resulting 
are given in Tables UTII and HVl respectively. 

Table III. The values of parameters in the best fit to the data up to the SPS energy (yT = 0.9TeV) 
in the analysis 2(fit to the data in lOGeV < k < Piarge = 4.3 x 10®GeV ). We obtain smaller 
error of ±*■"'"^(0) than in analysis 1 (Table HJ, since, as is seen in Eq. iim . ±(+>(0) has sizable 
effects only in the low energy region. For errors, see the caption in Table |I] 



C2 

Cl Co 

Pp' 

±(+>(0) 

Analysis 2 

0.0479 ± 0.0037 

-0.186 T 0.057 6.38 ±0.22 

7.26 T 0.33 

10.19 ± 1.72 


Table IV. The values of for the fit to data in lOGeV < k < Piarge = 4.3 x 10® GeV(Analysis 2). 
For Nf and Na{Np), see the caption in Table, im 



X^Nf xl/N^ 

Xl/Np 

Analysis 2 

14.1/37 8.8/32 

5.3/10 


The predictions at LHC energy {y/s = 14TeV) in terms of the best fit values of 
high-energy parameters in Table IThI are 


= 107.8 ± 2.4 mb, pPP = 0.129 ± 0.004 , 


(3-5) 












On the Discrepancy of pp, pp Total Cross Sections at = l.STeV 7 

where the errors correspond to the one standard deviation of C 2 . Essentially the same 
prediction are obtained as Eq. imi of the analysis 1, although the errors are slightly 
smaller. Our result is stable independently of the choices of the htting energy range. 

§4. Concluding remarks 

In §3, we have investigated a possibility to resolve the discrepancy between 
E710, E811 and CDF, using the experimental data of and up to the SPS 
experiments (y^ = O.DTeV). 

We came to the conclusion that only the data of E710 is consistent with the 
prediction, Eq. in the one standard deviation although we can exclude neither 
CDF nor E811 results in the two standard deviations. In our previous paper, ref. 1) 
we concluded that the precise measurements of in the coming LHC measurements 
will resolve this discrepancy at ^/s = l.STeV. It would still be worthwhile , however, 
to fix this problem in the CDF and DO experiments, since these values play an 
important role to search for and in the higher energy regions. 

Appendix A 

- Reanalysis of our predictions at the LHC (y/s=14TeV) with 

parameter - 

In our previous work,^) we exploited the experimental data and p^^^ above 

Piab=70GeV up to Tevatron energy {^/s = 1.8TeV) to predict and p^~^^ in the 

LHC region, based on Eq. (EIHI) of p^^\ not by Eq. (ESI). Although the effect of 
the parameter F^P{0) in the new formula (Eq. (12-91) 1 is not large in the high energy 
region, we show the results of the analyses based on Eq. ((231) here for completeness. 

Corresponding to ref. 1) two independent analyses are done: one includes the 
E710/E811 data at y^=1.8TeV denoted as fit 2 in ref. 1), and the other includes the 
CDF datum of at the same energy denoted as fit 3 in ref. 1). The results of the 
simultaneous fit to and p^^'^ are compared with the previous results^^ in Fig. ISl 
The fit to is slightly improved in the lower energy region, while the result of 
is almost the same as the previous one. The obtained values of parameters and 
the resulting are given in Table IVl and Table IVTl respectively. 

Table V. The best-fit values of parameters in the fit 2 (fit up to Tevatron-collider energy including 
E710/811 data) and fit 3 (including CDF datum). The errors here correspond to the one 
standard deviation of C 2 . (See the caption in Table 0.) 



C2 

Cl 

CO 

Pp' 

±(+>(0) 

fit 2 

0.0424 ± 0.0041 

-0.099 =F 0.069 

6.04 ±0.28 

7.61 ± 1.55 

12.48 ±0.73 

fit 3 

0.0496 ± 0.0043 

-0.205 =F 0.072 

6.44 ±0.29 

7.20 ±0.81 

12.78 ±0.72 


The fit to in the lower energy region is improved in comparison with the 
previous result, as can be seen in Fig. (21 Correspondingly much smaller is obtained 
in Table IVTl which is compared with the previous values, yp=8.4(6.9) for fit 2(3). 
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(a) : All region 


(b) 


A+). 


Low energy region 


mb 



-—^^^—' Ecm 

5 10 20 OOlOOTeV 



(c) : High energy region (d) 

Fig. 2. Predictions for and in terms of the Analysis 2. The fit is done for the data up to 
SPS energy, in the region 10GeV< k < 4.3x 10®GeV(4.54GeV < < 0.9TeV) which is shown by 

the arrow. For each figure, see the caption in Fig.0 The thin dot-dashed lines represent the one 
standard deviation of C 2 .(See the caption in TableQ]) The corresponding values of parameters are 
{c2,ci,co,(3p,,F<-+'>{0)) = (0.0479±0.0037,-0.186T0.056, 6.38T0.21, 7.26=f0.31, 10.19T0.31). 

Table VI. The values of for the fit 2 and fit 3. Np and Na-{Np) are the degree of freedom and 

the number of data points in the fitted energy region. 



xVNf 

Xl/N. 

Xl/Np 

fit 2 

11.6/22 

7.9/18 

3.7/9 

fit 3 

10.9/22 

8.7/18 

2.1/9 
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Predicted values of and at LHC energy(-y/i=14TeV) and at cosmic-ray 
energy {Piab=^ x 10^°eV) are given in Table IVTIl 

Table VII. The predictions of and at the LHC energy y/s = Ecm = 

14TeV(Pia6=1.04x 10®GeV), and at a very high energy Piab = 5 • 10^°eV (vT=i?cm=967TeV.) 
in the cosmic-ray region. The errors correspond to one standard deviation of C 2 . 



4J't^(vT=14TeV) 

p(+)(V7=14TeV) 

(rZ{Piab=5 ■ lO^OeV) 

P^+HTa6=5 ■ lO^OeV) 

fit 2 

104.2 ± 2.3mb 

0.123 ±0.004 

191 ± 8mb 

0.100 ± 0.003 

fit 3 

109.3 ± 2.4mb 

0.130 ±0.004 

206 ± 8mb 

0.105 ±0.003 


The predictions combining the two results in Table ivm are 

aZ = 106.8 ± S.lsyst ± 2.4tat mb, ffP = 0.127 ± O.OOTsyst ± 0.004stat 

(jfP = 198 ± 16syst ± 8stat mb, pPP = 0.103 ± 0.004syst ± O.OOS^tat (A-1) 

at the LHC energy(y^ = Ecm = 14TeV) and the cosmic-ray energy {Piab = 5 x 
lO^^eV), respectively. The above resnlts are almost the same as the previous ones, 
Eq. (13) of ref. 1). Here we obtain fairly large systematic uncertainty again coming 
from the data treatment at the Tevatron-energy. 
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mb 



(a) previous result^) 



mb 



(b) : using (0) parameter 



(c) previous result^) (d) using F(+^{0) parameter 

Fig. 3. Predictions for cr*-^^ and compared with the previous resnlts: The new results using 
^’^'*"^(0) parameter are shown by right figures, (b) and (d), respectively, which are compared 
with the left figures, (a) and (c), of the previous analyses. Predictions in terms of the fit 
2(3) are shown by green(blue) lines, and the thin dot-dashed lines represent the one standard 
deviation of C 2 . (See the caption in Table H]) The corresponding values of parameters are given 
in Table IVl 








































